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Abstract 

We develop the color dipole phenomenology of diffractive photo- and electropro- 
duction 7* N —> V(V') N of light vector mesons (V(1S) = (p°,uj ,p ) and their radial 
excitations (V'(2S) = <p' , u/, p'). The node of the radial wave function of the 2S states 
in conjunction with the energy dependence of the color dipole cross section is shown 
to lead to a strikingly different Q 2 and v dependence of diffractive production of the 
V(1S) and V'(2S) vector mesons. We discuss the restoration of flavor symmetry 
and universality properties of production of different vector mesons as a function of 
Q 2 + my. The color dipole model predictions for the p° and <jp production are in good 
agreement with the experimental data from the EMC, NMC, ZEUS and HI collabora- 
tions. We present the first direct evaluation of the dipole cross section from these data. 
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1 Introduction 



Diffractive electroproduction of vector mesons 

7 *p-y p , V = p°, uj°, 0°, J/ty, T (1) 

at high energy v offers a unique possibility of studying the pomeron exchange at high energies 
[H, ||, 0, f|, §, H 0, |] . Particularly important is the observation that the transverse size of 
the photon shrinks with the increase of its virtuality Q 2 . This property can conveniently be 
quantified in the mixed (r, z) lightcone technique [g, [IIJ, in which the high energy hadrons 
and photons are described as systems of color dipoles with the transverse size r frozen during 
the interaction process. Interaction of color dipoles with the target nucleon is quantified by 
the color dipole cross section o~{y,r) whose evolution with the energy v is described by the 
generalized BFKL equation [I^, |ll| (for a related approach see also |12| ). The shrinkage of 



the photon with Q 2 together with the small-size behavior of the dipole cross section (~ r 2 ) 
leads to what has come to be known as a scanning phenomenon [|l^, |4|, ||, ||: the V(1S) 
vector meson production amplitude is dominated by the contribution from the dipole cross 
section at the dipole size r ~ r^, where r$ is the scanning radius 

r s « . . (2) 

I m 2 v + Q 2 



This scanning property makes the vector meson production an ideal laboratory for testing 
the generalized BFKL dynamics [gBFKL hereafter]. At large Q 2 and/or for heavy vector 
mesons, the amplitude of reaction ([!]) becomes short- distance dominated and is perturba- 
tively calculable in terms of the short-distance behavior of the vector mesons wave function. 
However, the asymptotic short distance formulas 0, 0] are not yet applicable at the moder- 
ate Q 2 < 20 GeV 2 of interest in the present fixed target and HERA experiments where the 
scanning radius rs is still large due to a large scale parameter A w 6 in (fj) as derived in 0. 
For this reason, the onset of the short-distance dominance is very slow and there emerges a 
unique possibility of studying the transition between the soft and hard interaction regimes 
in a well controlled manner. Furthermore, the scanning phenomenon allows to directly test 
the steeper subasymptotic energy dependence of the dipole cross section at smaller dipole 
size r, which is one of interesting consequences of the color dipole gBFKL dynamics |Tl|, 



The scanning phenomenon has particularly interesting implications for the diffractive 
production of the 25 radially excited vector mesons 

^p^V'p, V'(2S)= P ',uj', V.. (3) 

Here one encounters the node effect: a tricky and strong cancellation between the large and 
the small size contributions to the production amplitude i.e., those above and below the node 
position r n in the 25 radial wave function 0, [13|, [15] respectively. The node effect is the 
only dynamical mechanism that gives a strong natural suppression of the photoproduction of 
excited vector mesons V'(25) vs. V(15) mesons. For instance, it correctly predicted 0, |T3| 
the strong suppression of real photoproduction of the ^' compared to the J/^f observed 
in the NMC experiment and confirmed recently in the high statistics E687 experiment 



| T7|j . In anticipation of the new experimental data on real and virtual V photoproduction 
from HERA, it is important to further explore the salient features of the node effect in the 
framework of the color dipole gBFKL dynamics. At moderate Q 2 , the scanning radius rs 
is comparable to r n . First, for this reason even a slight variation of rs with Q 2 leads to 
a strong change of the cancellation pattern in the V^'(25) production amplitude and to an 
anomalous Q 2 dependence for the electroproduction of the radially excited vector mesons 
P, [T3|, |T5|]. Second, the cancellation pattern is sensitive also to the dipole-size dependence 
of the color dipole cross section a(u, r) which in the gBFKL dynamics changes which energy 
v leading to an anomalous energy dependence for producing the V'(2S) vector mesons as 
compared to a smooth energy dependence for the ^(15) ground state vector mesons. This 
anomalous Q 2 and energy dependence of the V'(2S) production offers a unique signature 
of the 25 radial excitation vs. the D-wave state. Third, at very small Q 2 , the V'(2S) 
production amplitude can be of opposite sign with respect to that of the V^(15) production 
amplitude (the overcompensation scenario of ref. to then conform to the same sign at 
larger Q 2 (the undercompensation scenario of ref. ]15[). Here we wish to emphasize that 
the relative sign of the V and V production amplitudes is experimentally measurable using 
the so-called Soding-Pumplin effect ( |18|, [19| , see also []20f ) . 



In this paper we develop the color dipole phenomenology of diffractive photo- and elec- 
troproduction of the 15 ground state and of the 25 radially excited vector mesons. As stated 
above, for a large scanning radius, the large distance contribution to the production ampli- 



tude is not yet negligibly small in the so far experimentally studied region of Q 2 , in particular 
in the 2S meson production. In this paper we show that the Q 2 and energy dependence of 
the diffractive production of vector mesons offers a unique possibility of studying how the 
color dipole cross section changes from the large nonperturbative to the small perturbative 
dipole size. The problem can be attacked both ways. First, we present detailed predictions 
using the color dipole cross section ||, which gives a very good quantitative description 



of the proton structure function from very small to large Q 2 . Second, we can invert the 
problem and evaluate the color dipole cross section from the corresponding experimental 
data. Such an evaluation of the dipole cross section is presented here for the first time. 

The paper is organized as follows. In section 2 we formulate the color dipole factorization 
for vector meson production amplitudes. In section 3 we present our numerical results. We 
find good agreement with the experimental data from the fixed target and HERA collider 
experiments. The subject of section 4 is the anomalous Q 2 and energy dependence of 
electroproduction of 25* radially excited vector mesons. In section 5 we discuss the scaling 
relations between production cross sections for different vector mesons and the restoration 
of flavor symmetry in the variable Q 2 + rriy. We comment on how the scanning phenomenon 
enables a direct comparison of the spatial wave functions of the p° and a; mesons. The first 
evaluation of the dipole cross section from the experimental data is presented in section 6. 
In section 7 we summarize our main results and conclusions. In the Appendix we describe 
the lightcone parameterization of wave functions of V(1S) and V(2S) vector mesons used 
in our analysis. 

2 Color dipole factorization for vector meson produc- 
tion 

The Fock state expansion for the relativistic meson starts with the quark-antiquark state 
which can be considered as a color dipole. The relevant variables are the dipole moment 
r which is the transverse separation (with respect to the collision axis) of the quark and 
antiquark and z - the fraction of the lightcone momentum of the meson carried by a quark. 
The interaction of the relativistic color dipole with the target nucleon is described by the 
energy dependent color dipole cross section a(v, r). The many gluon contributions of higher 
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Fock states qqg... become very important at high energy v. The crucial point is that in the 
leading log- the effect of higher Fock states can be reabsorbed into the energy dependence 
of a{y, r), which satisfies the generalized BFKL equation ( |]T0| , |TT|1). The flavor blind (one 
should really say flavor tasteless) dipole cross section unifies the description of various 
diffractive processes. To apply the color dipole formalism to deep inelastic and quarkonium 
scattering and diffractive production of vector mesons one needs the probability amplitudes 
\l/ 7 * (f, z) and \lV (f, z) to find the color dipole of size r in the photon and quarkonium (vector 
meson), respectively. The color dipole distribution in (virtual) photons was derived in [p], |10[1 . 
In terms of these probability amplitudes, the imaginary part of the virtual photoproduction 
of vector mesons in the forward direction (t = 0) reads 

ImM = (V\a(v, r) | 7 *) = J dz J al 2 r a(u, r)^* v (r, z)^ r (r, z) (4) 

whose normalization is do~/dt\ t=0 = \A4\ 2 /16ir. For small size heavy quarkonium the prob- 
ability amplitude \lv(r,,2) can safely be identified with the constituent quark-antiquark 
quarkonium wave function. The color dipole factorization (|J) takes advantage of the diago- 
nalization of the scattering matrix in the (r, z) representation, which clearly holds even when 
the dipole size r is large, i.e. beyond the perturbative region of short distances. Due to this 
property and to the fact that in leading log^ the effect of higher Fock states is reabsorbed in 
the energy dependence of the dipole cross section cr(z/, r), as a starting approximation we can 
identify the probability amplitude \Ev(r*, z) for large size dipoles in light vector mesons with 
the constituent quark wave function of the meson. This provides a viable phenomenology of 
diffractive scattering which is purely perturbative for small size mesons and/or large Q 2 and 
small scanning radius r$ and allows a sensible interpolation between soft interactions for 
large dipoles and hard perturbative interactions of small dipoles. For light quarkonia and 
small Q 2 , this implies the assumption that small-size constituent quarks are the relevant 
degrees of freedom and the spatial separation of constituent quarks is a major dynamical 
variable in the scattering process. f\ The large-r contribution to the production amplitude 
(|^) depends on the both dipole cross section for large-size dipoles and the amplitudes of 

distribution of large-size color dipoles and/or the nonperturbative wave functions of light 

* See also earlier works on the color dipole analysis of hadronic diffractive interactions which used 
constituent quark wave functions for the color dipole distribution amplitudes [^2| . 
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vector mesons at large r, both of which are poorly known at the moment. Still, testing 
the predictions from such a minimal model is interesting for its own sake and can shed a 
light on the transition between the soft and hard scattering regimes which is still far from 
understood. An analysis of sensitivity to models of the nonperturbative wave functions of 
vector mesons and of how one can disentangle the effects of large r behavior of the wave 
function and of the dipole cross section, goes beyond the scope of the present exploratory 
study. 

The energy dependence of the dipole cross section is quantified in terms of the dimen- 
sionless rapidity £, which in deep inelastic scattering equals £ = log^. Considerations of 
intermediate masses in diagrams for exclusive production of vector mesons show that to the 
considered leading log- approximation one must take £ = log^— , where 

X x ef f 

Q 2 + m 2 v 

and my is a mass of the vector meson. The pomeron exchange dominance holds when the 
Regge parameter is large, 

1 2m v v „ 

— = iro ■ » 1 - ( 6 ) 



x e ff (Q 2 + m 



Hereafter we write the amplitudes in terms of a(x e ff, r). The spin independence of the dipole 
cross section a(x e ff,r) in (|J) leads to the s-channel helicity conservation: the transversely 
polarized photons produce transversely polarized vector mesons and the longitudinally po- 
larized vector mesons are produced by longitudinal (to be more precise, the scalar one) 
photons. More explicitly, the form of the forward production amplitudes for the trans- 
versely (T) and the longitudinally (L) polarized vector mesons in terms of the lightcone 
radial wave function 0(r, z) of the qq Fock state of the vector meson reads || 

lmM T {x eff ,Q ) = - 



(27T) 2 

d 2 va(x eff ,r) * - {m 2 q K (er)(j)(r, z) - [z 2 + (1 - z) 2 )eKx{er)d r 4){r, z)} 

dr 2 <j(x. 



(m 2 v + Q 2 ) 2 

N /Zttv 

lmM L (x eff ,Q 2 ) 



J d y-^W T (Q 2 ,r 2 ) (7) 



(2tt) 2 my 
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d 2 ra(x e ff,r) J dz {[m 2 , + z(l — z)rriy]Ko(er)(f)(r, z) — eKi(er)d r (j)(r, z)} 

= I h^L f ^1^uA Wl (Q 2 r 2 ) (8) 

(m^ + Q 2 ) 2 m v J r 2 r 2 lW ' ] K ) 

where 

e 2 = m\ + z{l - z)Q 2 , (9) 

a em is the fine structure constant, N c = 3 is the number of colors, Cy = §, § for 

p°, u; , 0°, J/\l/ production, respectively and Kqj_(x) are the modified Bessel functions. The 
detailed discussion and parameterization of 0(r, z) is given in the Appendix, here we only 
mention that the form of 0(r, z) we use has the hard-QCD driven short distance behavior 
and gives the electromagnetic form factor of mesons which has the correct QCD asymptotic 
behavior. At large r we follow the conventional spectroscopic models []23] and constrain 
the parameters of the wave functions by the widths of the leptonic decays V,V — ► e + e~, 
the radii of the vector mesons and the 2S-1S mass splitting. The terms oc 0(r, z)Ko(er) 
and oc d r (j)(r, z)eK 1 (er) 1 i.e., d r (j)(r, z)d r K (er), in the integrands of (0) and (§) derive from 
the helicity conserving and helicity nonconserving transitions 7* — > qq and V — > qq in the 
A^-f^ and V^^^ vertices (see Bjorken et al. [^4]] and 0, the technique of calculation 



of traces in the spinorial representation of the relevant Feynman amplitudes is given in 
II and need not be repeated here; for the related Melosh transformation analysis see the 



recent Ref. [25]). The latter are the relativistic corrections, for the heavy quarkonium 
the nonrelativistic approximation has a rather high accuracy, the relativistic corrections 
become important only at large Q 2 and for the production of light vector mesons. Eqs. 

give the imaginary part of the production amplitudes; one can easily include small 
corrections for the real part by the substitution p6|, 



a(x e// , r) ( 1 - i ■ J • ^-j— ] cr(x eff , r) 

\ 2 d\ogx eff ) 



l-i- a v (x eff ,r) 



<?(xeff,r) (10) 



For brevity, in the subsequent discussion we suppress the real part of the production ampli- 
tude; it is consistently included in all the numerical calculations. 

The color dipole cross section is flavor blind. The only kinematical sensitivity to the 
vector meson produced comes via the rapidity variable, see Eq. @. For small r, in the 
leading log- and leading log^-, i.e., leading logQ 2 , the dipole cross section can be related 



| 27fl to the gluon structure function G(x, Q 2 ) of the target nucleon through 



a(x,r) = ^-r 2 a s (r)G(x,Q 2 ), (11) 

where the gluon structure function enters at the factorization scale Q 2 ~ ^ (for the origin 
of the large scale factor B ~ 10, see [28]). The integrands of (0),(||) are smooth at small 
r and decrease exponentially at r > 1/e due to the exponential decrease of the modified 
Bessel functions. Together with the oc r 2 behavior of the color dipole cross section (|TT1), 
this implies that the amplitudes (0),(§) receive their dominant contribution from r « rs- 
(Eq. (|2]) assumes that the scanning radius rs is substantially smaller than the radius Ry of 
the vector meson.) Then, a simple evaluation gives f| 

lmM T oc r 2 s a(x eff , r s ) oc a{x eff , r s ) oc (g2 + m ^ 2 ( 12 ) 

and 



lmML " 1^ Mt x ^ a(x '"- rs) K l!^ (Qi TK? (13) 

respectively. f\ The prediction of the dominance of the longitudinal cross section at large Q 2 
is shared by all the models of diffractive leptoproduction, starting with the vector dominance 
model (H, [5], [/J, for the excellent review of early works on photo- and electroproduction of 
vector mesons and on vector dominance model see Bauer et al. PD| ) and is confirmed by 
all the experiments on leptoproduction of the p° at large Q 2 |30], |3l], [32| . The first factor 
oc r| oc 1/(Q 2 +rriy) in (|T^) comes from the overlap of wave function of the shrinking photon 
and that of the vector meson. The familiar vector dominance model (VDM) prediction is 
Mf oc ^ m 2 1 + Q2} &toti.pN) , whereas in our QCD approach a small u(x e ff,r s ) oc r| oc l/(Q 2 + 
y) enters instead of a to t(yN). In (p^2|),(p^3|) we show only the leading Q 2 dependence, 



in 



suppressing the phenomenologically important departure form the law a(x, r) oc r 2 , whose 
large Q 2 dependence can be related to scaling violations in the gluon density (see ( |TTD and the 
discussion below). We recall that the shrinkage of the virtual photons and/or the decrease 
of the scanning radius rs with Q 2 is the origin of color transparency effects in diffractive 

leptoproduction of vector mesons off nuclei 0, f|, [5], |33] . The important confirmation of the 

' Unless otherwise specified, for each flavor, my will always be the mass of the ground state IS vector 
meson. 
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quantitative predictions ||, |[ of color transparency effects based on the same technique as 



used here came from the E665 experiment p4 



More accurate analysis of the scanning phenomenon can be performed in terms of the 



weight functions Wt,l(Q 2 , f 2 ) which are sharply peaked at r « At,l/ \jQ 2 + m \i i n ^ ne rele- 
vant variable logr the width of the peak in Wl{Q 2 , r 2 ) at half maximum equals A logr « 1.2 
for the J/^f production and A logr « 1.3 for the p° production and varies little with Q 2 ||. 
The values of the scale parameter Atl turn out to be close to A ~ 6, which follows from 
r s — 3/e with the nonrelativistic choice z = 0.5; in general At,l > 6 and increases slowly 
with Q 2 . This Q 2 dependence of At,l comes from the large-size asymmetric qq configurations 
when, for instance, the antiquark and the quark in the photon and in the vector mesons carry 
a very large and a very small fraction of the meson momentum respectively (or the other way 
around). A comparison of the integrands in eqs. (|7|) and (^j) shows that the latter contains 
an extra factor z(l — z) which makes considerably smaller the contribution from asymmet- 
ric configurations to the longitudinal meson production. For completeness, we quote the 
results of §: A T , L {J/^]Q 2 = 0) « 6, A T:L (J/^;Q 2 = 100 GeV 2 ) « 7, A L (p°;Q 2 = 0) « 
6.5, A L {p°;Q 2 = 100 GeV 2 ) w 10, A T {p°;Q 2 = 0) w 7, A r (p°,Q 2 = 100 GeV 2 ) w 12. 

An alternative formulation of the slow onset of the purely perturbative regime can be 
seen as follows: at very large Q 2 when the scanning radius is very small, the dipole cross 
section a(x e ff, r) and the vector meson production amplitudes are proportional to the gluon 
density G(x e ff, Q 2 ) at the factorization scale Q 2 = t(Q 2 + my) (see also Refs. 0, 0] which 
use a different technique of the momentum-space wave functions, related to the color dipole 
factorization by the Fourier-Bessel transform; the detailed comparison with the work of 
Brodsky et al. will be presented below in Section 6). The large values of At,l previously 
quoted, reflect into very small values of r [H: in the interesting region of Q 2 > 10 GeV 2 
one finds t t ,l{J/^) ~ 0.2, r L (p°) w 0.15 and r T (p°) ^0.07-0.1, which is different and 
substantially smaller than the values r = 0.25 suggested in and r = 1 suggested in 
0. Very large Q 2 values are needed for reaching the perturbatively large Q 2 and for the 
applicability of the pQCD relationship ([□])• 

Consequently, for the domain presently under experimental study, Q 2 +my <10-20 GeV 2 , 
the production amplitudes receive substantial contribution from semiperturbative and non- 
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perturbative r. In [[H], this contribution was modeled by the energy independent soft 
cross section a^ npt \r). The particular form of this cross section successfully predicted [^T] 
the proton structure function at very small Q 2 recently measured by the E665 collaboration 



p5| and also gave a good description of real photoabsorption ||. As an example, in Fig. 1 



we present an evaluation of the vector meson-nucleon total cross section 

vtottYN) = ^ [ z2{ * Z _ z)2 J d 2 r {rnfrir, z) 2 + [z 2 + (1 - zf] [d r <J>(r, z)] 2 } a(x eff , r) .(14) 

The total cross section a tot (p°N) so found, is close to cr tot (7iN), and the rise of a tot (VN) 
with the c.m.s energy W is consistent with the observed trend of the hadronic total cross 
sections |3f|. In the color dipole picture the smaller values of a to t((pN) and <J tot ((f)'N) derive 



from the smaller radius of the ss quarkonium. In the simple model pi], || the rise of 
VtotiVN) is entirely due to the gBFKL rise of the perturbative component a^ pt \x e ff^r) of 
the dipole cross section. The rate of rise is small for two reasons: i) at moderate energy, 
a^ ptS) (x e ff,r) at large r is much smaller than the soft cross section a^ npt \r), ii) at large r 
the subasymptotic effective intercept of the gBFKL pomeron is small |TT], [TJ|]. The detailed 
description of the dipole cross section used in the present analysis is given in |2l|, |6| and 
will not be repeated here. It is partly shown below in Fig. 16. The reason why we focus 
here on this particular model is that its success in phenomenological applications makes 
it a realistic tool for the interpolation between soft and hard scattering regions. Once the 
vector mesons wave functions are fixed from their spectroscopic and decay properties, all 
the results for diffractive real and virtual photoproduction of vector mesons to be reported 
here do not contain any adjustable parameters. 

3 Diffractive p° and (f)° production: predictions and 
comparison with experiment 



The most interesting prediction from the color dipole dynamics is a rapid decrease of produc- 
tion amplitudes (JT^),(JT^) at large Q 2 . The broadest region of Q 2 was covered in the recent 



NMC experiment |3l] where special care was taken to minimize the inelastic production 
background which plagued earlier data on p° and J/\& production. In Fig. 2 we compare 
our predictions for p° and 0° production with the NMC data and the data from the HERA 
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experiments [p2| , |37| . Shown here is the observed polarization- unseparated cross section 
a(7* V) = 0^(7* — > V) + ecr L (7* — ► V) for the value of the longitudinal polarization e 
of the virtual photon taken from the corresponding experimental paper (typically, e ~ 1). 
The quantity which is best predicted theoretically is da(Y — > V)/dt\ t=0 ] in our evaluations 
of the total production cross section 17(7* — ► V) = -8(7* — > V)da(j* — > V)/dt\ t =o we use 
the diffraction slope -8(7* — > V) given in the corresponding experimental paper. 

Eqs. (0),© describe the pure pomeron exchange contribution to the production ampli- 
tude. While at HERA energies secondary Reggeon exchanges can be neglected since the 
Regge parameter uj is a very large, at the lower energy of the NMC experiment, (y) =(90- 
140) GeV, the Regge parameter uj is small and non-vacuum Reggeon exchange cannot be 
neglected. The fit to (?tot{lP) can, for instance, be cast in the form 

cr to t(jp) = &tp(ip) ■ + -jx) ( 15 ) 
where the term A/uj a in the factor / = 1 + A/uj a represents the non- vacuum Reggeon 
exchange contribution. The Donnachie-Landshoff fit gives A = 2.332 and A = 0.533 [pEfl . 
We do not know how large this non- vacuum contribution to p° production is at large Q 2 ; 
for a crude estimation we assume the Reggeon/pomeron ratio to scale with uj, which is not 
inconsistent with the known decomposition of the proton structure function into the valence 
(non-vacuum Reggeon) and sea (pomeron) contributions. Then, for the NMC kinematics 
we find / = 1.25 at uj ~ 70, Q 2 = 3 GeV 2 and / = 1.8 at uj ~ 9,Q 2 = 20 GeV 2 . This 
departure of / from unity provides a conservative scale for the theoretical uncertainties at 
moderate values of uj. Anyway, the Q 2 dependence of the Reggeon correction factor / is 
weak compared with the very rapid variations of M.t and M.l with Q 2 . The correction 
for the secondary exchanges, 7(7* — > p°) = / 2 crip(7* — * p°), brings the theory to a better 
agreement with the NMC data. The dipole cross section of 0, U correctly describes the 
variation of the p° production cross section by 3 orders in magnitude from Q 2 = to 
Q 2 = 16.5 GeV 2 . For 0° production, / = 1 due to the Zweig rule and the pure pomeron 
contribution correctly reproduces the magnitude of 17(7* — > 0°) and its variation by nearly 
three orders in the magnitude from Q 2 = to Q 2 = 11.3 GeV 2 . 

The specific prediction from the gBFKL dynamics is a steeper subasymptotic growth 
with energy of the dipole cross section a{y,r) at smaller dipole size r, which by virtue of 
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the scanning phenomenon translates into a steeper rise of (7(7* — > V) at higher Q 2 and/or 
for heavy quarkonia. This consequence of the color dipole dynamics was first explored in 
H; the p° wave function parameters used in || are slightly different from those used here 
but the difference in cr(7* — > p°) is marginal. The agreement of our high-energy results with 
the HERA data is good for both Q 2 = (Fig. 3) and large Q 2 (Fig. 2) and confirms the 
growth of the dipole cross section with energy expected from the gBFKL dynamics. 

The above high-Q 2 data are dominated by the longitudinal cross section; real photo- 
production (Q 2 = 0) measures the purely transverse cross section. In Fig. 3 we present 
our results with and without secondary Reggeon corrections (dajp^ — > p°)/dt\ t=Q and 
rfcr(7 — > p°)/dt\ t=0 = f 2 dujp(j — > p°)/dt\ t =o respectively) as a function of energy. The 
Reggeon correction factor f 2 brings the theory to a better agreement with the low energy p° 
production data ]38[ . Real photoproduction of p° is dominated by the soft contribution, the 



growth of the production cross section is driven by the rising gBFKL component of the dipole 
cross section. Our predictions for high energy agree well with the recent ZEUS data |39|. [HJ. 
The 0° production is pomeron dominated which implies / = 1. We find good agreement 
with the fixed target |4l| and ZEUS (|2 data on real photoproduction of the 0°, although 



the error bars are large (Fig. 4). Because in 0° photoproduction the relevant dipole sizes 
are smaller than in the p° case, (see the radii of p° and 0° in Table 1), we predict a steep en- 
ergy dependence of the 0° production forward cross section: da{^f — > 4>°)/dt\ t =o is predicted 
to grow by a factor rs 2.5 from 3.75 pb/GeV 2 at v = 175GeV, i.e., W = 18GeU, up to 
~ 8.84 pb/GeV 2 at W = 170 GeV at HERA. At W = 70 GeV we have a( 7 -> 0°) = 0.87/1& 
which agrees with the first ZEUS measurement (7(7 — ► 0°) = 0.95 ± 0.33 pb |42|]. More 
detailed predictions for the energy and Q 2 dependence of da{j* — > V)/dt\ t =o are presented 
in Fig. 5 and clearly show a steeper rise with energy at larger Q 2 (see also 0). 
In Fig. 6 we show our predictions for 

RLT ~ Q 2 - da T (r ^ V) ■ (16) 
The steady decrease of Rlt with Q 2 which implies a diminution of the dominance of the 

longitudinal cross section is a very specific prediction of the color dipole approach. It 

follows from a larger contribution from large size dipoles to the production amplitude for 

the transversely polarized vector mesons and larger value of the average scanning radius, i.e., 
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At ^ Al 0- This prediction can be checked with the higher precision data from HERA; 



the available experimental data [|3(], |3T], ^] agree with Rlt < 1 but have still large error 
bars. 

The Q 2 dependence of the observed polarization-unseparated cross section depends on 
the longitudinal polarization e of the virtual photon. To a crude approximation the color 
dipole dynamics predicts 

a(;' —*V)= , T (Y -+ V) + -+ V) cc ^,1^ (l + (17) 

If one approximates (|TT| ) by the (Q 2 + m v )~ n behavior, one finds n ~ 3 vs. n ~ 1 in the 
naive VDM. In (|l7|) we suppressed the extra Q 2 dependence which at large Q 2 comes from 
the scaling violations in the gluon density factor oc G 2 (x,r(Q 2 + my)), see ([TT|). For these 
scaling violations, at fixed x e ff and asymptotically large Q 2 we expect n < 3. In Fig. 7a we 
present our predictions for p° and 0° production at W = 100 GeV as a function of Q 2 + 



assuming for the longitudinal polarization e = 1 as in the ZEUS kinematics |32|. These 
cross sections can be roughly approximated by the oc (Q 2 + m\r)~ n law with the exponent 
n 2.4 for the semiperturbative r$ region 1 < Q 2 < 10 GeV 2 . At fixed W, x e ff varies 
with Q 2 and for the x e ff dependence of a(x e ff, Q 2 ) we predict n w 3.2 for the perturbative 
15 < Q 2 < 100 GeV 2 where rs is small. We strongly urge a careful analysis of the Q 2 
dependence in terms of the natural variable Q 2 + my (for more discussion see section 5 
below). For the sake of completeness, in Fig.5 we present also our predictions for the energy 
dependence of the polarization-unseparated production cross section a = a? + ea^ for the 
typical e = 1. 

4 Anomalies in electroproduction of 2S radially ex- 
cited vector mesons 

Here the keyword is the node effect - the Q 2 and energy dependent cancellations from 
the soft (large size) and hard (small size) contributions to the production amplitude of 
the V'(2S) radially excited vector mesons. When the value of the scanning radius rs is 
close to the node r n ~ Ry, these cancellations must exhibit a strong dependence on both 
Q 2 and energy due to the different energy dependence of the dipole cross section at small 
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(r < Ry) and large (r > Ry) dipole sizes. It must be made clear from the very beginning 
that when strong cancellations of the large and small region contributions are involved, 
the predictive power becomes very weak and the results strongly model dependent. Our 
predictions for the production of the V'(2S) radial excitations which we report here serve 
mostly as an illustration of the unusual Q 2 and energy dependence possible in these reactions. 
(Manifestations of the node effect in electroproduction on nuclei were discussed earlier, see 
H and @) 

In the nonrelativistic limit of heavy quarkonia, the node effect will not depend on the 
polarization of the virtual photon and of the produced vector meson. Not so for light vector 
mesons. The wave functions of the transversely and longitudinally polarized photons are 
different, the regions of z which contribute to the M.t and M.l are different, and the Q 2 
and energy dependence of the node effect in production of the transverse and longitudinally 
polarized V'(2S) vector mesons will be different. 



Let us start with the transverse amplitude. Two cases can occur ||I5| , the undercompen- 
sation and the overcompensation scenario. In the undercompensation case, the production 
amplitude (2S\a(x e ff,r)\j*) is dominated by the positive contribution coming from r <r n 
and the V(1S) and V'(2S) photoproduction amplitudes have the same sign. With our model 
wave functions this scenario is realized for transversely polarized p' and 0' (we can not, 
however, exclude the overcompensation scenario). As discussed in ||15|| , in the undercom- 
pensation scenario a decrease of of the scanning radius with Q 2 leads to a rapid decrease of 
the negative contribution coming from large r > r n and to a rapid rise of the V'(2S) /V(1S) 
production ratio with Q 2 . The stronger the suppression of the real photoproduction of the 
V'{2S) state, the steeper the Q 2 dependence of the V'{2S) /V{1S) production ratio expected 
at small Q 2 . With our model wave functions, the p'(2S)/p° and <p'(2S)/ ' (jP production ratios 
for the transverse polarization are predicted to rise by more than one order of magnitude 
in the range Q 2 < 0.5 GeV 2 , see Fig. 8; the V(2S) and V(1S) production cross sections 
become comparable at Q 2 > lGeV 2 , when the production amplitudes are dominated by 
dipole size r <r n |L5|, ||. 

For the longitudinally polarized p'{2S) and <p'{2S) mesons, our model wave functions 
predict overcompensation; at Q 2 = OGeV 2 the amplitude is dominated by the negative 
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contribution from r > r n . Consequently, with the increase of Q 2 , i.e. with the decrease 
of the scanning radius rs, one encounters the exact cancellation of the large and small 
distance contributions. Our model wave functions lead to this exact node effect in the 
dominant imaginary part of the production amplitude at some value Q 2 n ~ 0.5 GeV 2 for 
both the p' L (2S) and 4>' L {2S) production (see Fig. 6). The value of Q 2 n is slightly different for 
the imaginary and the real part of the production amplitude but the real part is typically 
very small and this difference will be hard to observe experimentally. Here we can not insist 
on the precise value of Q 2 n which is subject to the soft-hard cancellations, our emphasis is 
on the likely scenario with the exact node effect at a finite Q 2 n . 

We wish to emphasize that only the experiment will be able to decide between the 
overcompensation and undercompensation scenarios. For instance, let the p° and p'(2S) 
be observed in the tttt photoproduction channel. The Soding-Pumplin effect of interference 
between the direct, non-resonant 7p — > imp production and the resonant 7p — > p°(p')p — > 
wp production amplitudes leads to the skewed p° and p' mass spectrum. The asymmetry of 
the p°(p') mass spectrum depends on the sign of the p°(p') production amplitudes (|fL&[|, the 



detailed theory has been worked out in |19|]). The Soding-Pumplin technique has already 
been applied to the //(1600) mass region in jp — > 7r + 7r~ at 20 GeV studied in the SLAC 
experiment pDJ. Their fit to the p'(1600) mass spectrum requires that the sign of the p' 
production amplitude be negative relative to that of the p. Although the interpretation of 
this result is not clear at the moment, because there are two p'(1450) and p'(1700) states 
which were not resolved in this experiment, the Soding-Pumplin technique seems promising. 

With the further increase of Q 2 and decrease of the scanning radius one enters the 
above described undercompensation scenario. Although the radii of the ss and uu, dd vector 
mesons are different, the Q 2 dependence of p'(2S) / p° and <f>'(2S) / (fP production cross section 
ratios will exhibit a similar pattern. For both the transverse and longitudinally polarized 
photons, these ratios rise steeply with Q 2 on the scale Q 2 ~ 0.5 GeV 2 . At large Q 2 where the 
production of longitudinally polarized mesons dominates, the p'(2S)/p° and (f>'(2S)/4>° cross 
section ratios level off at ~ 0.3 (see Fig. 8). This large-Q 2 limiting value of the p'(2S)/p° and 
4>'{2S) 1 4>° cross section ratios depend on the ratio of V'(2S) and V(1S) wave functions at the 
origin, which in potential models is subject to the detailed form of the confining potential 
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| 23| . It is interesting that due to the different node effect for the T and L polarizations, we 
find R LT (2S) < Rlt(1S) , see Fig. 6. 

In Fig. 7b we present our predictions for the Q 2 dependence of the polarization-unseparated 
cross section (7(7* -> V'(2S)) = <x T (7* -> V'(2S)) + ea L (Y -> 1/'( 2 ^)) at the HERA 
energy = 100 GeV^ assuming e = 1. In Fig. 9 we show the Q 2 dependence of the 
polarization-unseparated forward cross section ratios da{Y — > p'(2S))/da( , ~f* — > p°) and 
dcr(7* — > (j)'(2S))/da(j* — >■ 0°) at = 100 GeV. Due to its smallness, the anomalous prop- 
erties of &l(2S) at small Q 2 are essentially invisible in the polarization-unseparated V'(2S) 
production cross section shown in Figs. 7b, 9 and 10. In contrast to (1(7* — > V(1S)), which 
falls monotonically and steeply from Q 2 = OGeV 2 on, the £"(7* — > V'(2S)) shown in Fig. 7b 
exhibits a weak rise at small Q 2 . At Q 2 large enough that the scanning radius r$ < Rv and 
the node effect becomes negligible, we predict very similar dependence on Q 2 + rriy of the 
V'(2S) and V(1S) production. 

Color dipole dynamics uniquely is the source of such a tricky Q 2 dependence of the 
V'(2S)/V(1S) production ratio. We already mentioned about the experimental confirmation 
|T6| . [17| of the node effect predicted in ^' production ||. Further experimental confirmations 
of the node effect, in particular of the unique overcompensation scenario which is possible 
for light vector mesons, would be extremely interesting. The available experimental data on 
real photoproduction of radially excited light V'{2S) mesons confirm (7(7 — > V {2S)) / a{j — > 
V(1S)) 1, but are still of a poor quality and the branching ratios of the V decays are 
not yet established (for the review see [H3fl and the Review of Particle Properties [fl5|]). For 



instance, the FNAL E401 experiment at v « 100 GeV found (7(7 -> ^'(1700, = 
8.0 ± 2.7(stat) ± 1.4(s?/st) nb to be compared with 17(7 — > (/)) « 0.55 /ib (see Fig. 4). In the 
p family, the very spectroscopy of the p' mesons is not yet conclusive ffl, There are 
two p' states, p'(1450) and p'(1600), the 2S and D-wave assignment for these states is not 
yet clear. The first high energy data on the p'(1450) and p'(1700) leptoproduction were 
reported by the E665 collaboration These E665 data refer to the coherent production 
on Ca target. For the p'(1700), they exhibit a strong rise of R21 = u{p' — > 47r)/cr(p — > 2tt) 
with Q 2 by more than one order in magnitude from (0.004 ± 0.004) at Q 2 = 0.15 GeV 2 
to (0.15 ± 0.07) at Q 2 = 4.5 GeV 2 . Such a steep Q 2 dependence is perfectly consistent 
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with our expectations for the production of radially excited 2S light vector mesons. For the 
p'(1450) there is a weak evidence of a nonmonotonic Q 2 dependence: R 2 \ = (0.035 ± 0.011) 
at Q 2 = 0.15 GeV 2 followed by decrease down to R 2 i = (0.012 ± 0.004) at Q 2 = 0.3 GeV 2 
and then to an increase and leveling off to R 2 i = (0.08 ± 0.04) at larger Q 2 > 2 GeV 2 . Such 
a Q 2 dependence of R21 would be natural for a D-wave state which has a nodeless radial 
wave function. If these E665 observations will be confirmed in higher statistics experiments, 
then the color dipole interpretation of the Q 2 dependence would strongly suggest the 2S and 
D-wave state assignments for the p'(1700) and p'(1450), respectively. We remind the reader 
that, for a quantitative comparison with the predictions of our model (shown in Fig. 10), 
the E665 results for R 2 ± must be corrected for the branching ratio B(p' — > Air), which is still 



experimentally unknown [ [|5] . 

The energy dependence of the p'{2S), 4>'{2S) real photoproduction is shown in Fig. 10 and 
has its own peculiarities. In the color dipole gBFKL dynamics, the negative contribution 
to the 25* production amplitude coming from large size dipoles, r > r n , has a slower growth 
with energy than the positive contribution coming from the small size dipoles, r > r n . 
For this reason, in the undercompensation regime the destructive interference of the two 
contributions becomes weaker at higher energy and we predict a growth of the V'(2S) /V(1S) 
cross section ratios with energy. Taking only the pure pomeron contributions into account, 
we find for the forward cross section ratio da{^ p'(2S)/da(>y p°) = 0.041 at W = 
15 GeV, which at HERA energies increases to 0.063 and 0.071 at W = 100 GeV and W = 
150 GeV, respectively. Whereas in p and p' production one must be aware of the non- 
vacuum Reggeon exchange contributions at lower energy, in the pomeron dominated <//, 4> 
real photoproduction we find a somewhat faster rise of da^ — > <f)'(2S)/da('y — > 0°) with 
energy from 0.054 at W = 15 GeV to 0.089 and 0.099 at W = 100 GeV and W = 150 GeV, 
respectively. 

If the leptoproduction of the longitudinally polarized V^S 1 ) will be separated experi- 
mentally, we will have a chance of studying the Q 2 and energy dependence in the overcom- 
pensation scenario. Start with the moderate energy and consider Q 2 very close to Q 2 n but 
still < Q 2 n . In this case the negative contribution from r > r n takes over in the V[(2S) 
production amplitude. With increasing energy, the positive contribution to the produc- 
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tion amplitude rises faster and ultimately takes over. At some intermediate energy, there 
will be an exact cancellation of the two contributions to the production amplitude and the 
longitudinal V[(2S) production cross section shall exhibit a minimum at this energy (the 
minimum will partly be filled because cancellations in the real and imaginary part of the 
production amplitude are not simultaneous). With our model wave functions, we find such 
a nonmonotonic energy dependence of the p' L (2S) and (f>' L (2S) production at Q 2 ~ 0.5 GeV 2 , 
which is shown in Figs. 8 and 10. At higher Q 2 and smaller scanning radii r$ the energy 
dependence of V[(2S)/Vl(1S) production ratio becomes very weak. 

Finally, a brief comment on the t-dependence of the differential cross sections is in or- 
der. For the IS vector mesons we expect the conventional diffractive peak with smooth and 
gentle energy dependence. For the radially excited vector mesons the t-dependence can be 
anomalous. The point is that the large size contribution to the V'(2S) meson production 
amplitude has steeper ^-dependence that the small size contribution. The destructive in- 
terference of these two amplitudes can lead to two effects: i) the diffraction slope in the 
V'(2S) meson production will be smaller than in the V(1S) meson production, ii) the effec- 
tive diffraction slope for the V'(2S) meson production decreases towards small t contrary 
to the familiar increase for the V(1S) meson production. High statistics data on the p', <ft' 
production at HERA are needed to test these predictions. More detailed discussion of the 
diffraction slope will be presented elsewhere. 

5 Scaling relations between production of different vec- 
tor mesons 

The color dipole cross section is flavor blind and only depends on the dipole size. The results 
(|l2l) , (PH) for the production amplitudes strongly suggest the restoration of flavor symmetry, 
i.e., a similarity between the production of different vector mesons when compared at the 
same value of the scanning radius rs and/or the same value of Q 2 +my. f\ Such a comparison 
must be performed at the same energy, which also provides the equality of x e ff at equal 

Q 2 + m v . Evidently, the value of Q 2 must be large enough so that the scanning radius rs 

tFor the first considerations of the restoration of flavor symmetry in diffractive production of vector 
mesons off nuclei see || , the scaling relations between diffraction slopes for j*p — ► Vp are discussed in M . 
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is smaller than the radii of vector mesons compared. 

In order to illustrate the above point we present in Figs. 11 and 12 the ratio of 
forward production cross sections R((J/^/)/p°;Q 2 ) = da(j* — > J/^)/da(j* — > p°) and 
R(<f)°/p°; Q 2 ) = der(7* — > 0°) / 'da{pf* —> p°) as a function of the c.m.s energy W at different 
Q 2 (here we use for the J/ty production cross section the values obtained from a recent cal- 



culation ||48|| , which practically coincide with those of ref.|6[], the slight difference being due 
to a somewhat different J/^> wave function). Here we compare the polarization- unseparated 
cross sections a = o T + tcr L , taking for the definiteness e = 1 which is typical of the HERA 
kinematics. These ratios exhibit quite a strong Q 2 dependence, which predominantly comes 
from the Q 2 dependence of the factor 



' Q 2 + m 2 



Vi 



Q 2 + m 2 



"v 2 / 

which changes rapidly when the two vector mesons have different masses. The energy 
dependence of the cross section ratios taken at the same Q 2 derives from the different 



energy dependence of the dipole cross section which enters at different radii r$i ■ ~ r — 

y/Q 2 

in the numerator and denominator of the V\jV% cross section ratio, 



a(7* -> V 2 ) a 2 (u, r S2 ) 
In the HERA energy range we predict R((J/^)/p°;Q 2 = 0) = (7(7 -> J/^)/a{^ -> p°) = 
0.022 at W = 70GeV and 0.028 at W = 150 GeV, which agrees with the experimentally 
observed ratio 0.0034 ± 0.0014 of HI (W ~ 70 GeV) |3J and 0.0045 ± 0.0023 of ZEUS 
(W = 150 GeV) [[H H, 0. Notice the rise of R((J/V)/p°; Q 2 ) by more than 3 orders in 
the mag nitude from Q 2 = to Q 2 = 100 GeV 2 . Our result for the ratio R((f)°/p°;Q 2 = 
0) = da{pf — > 0°)/c?<Tip(7 ~~ * P°) shown in Fig. 11 is substantially smaller than the factor 2/9 
expected from the naive VDM, in a very good agreement with the experiment ( |[41|| and 
references therein). This suppression is a natural consequence of the color dipole approach 
and derives from the smaller radius of the ss quarkonium and smaller transverse size of the ss 
Fock state of the photon as compared to the radius of the p° and size of the uu, dd Fock states 
of the photon, respectively, cf. Table 1. For increasing Q 2 s, the ratio R{<p / p°; Q 2 ) overshoots 
the VDM ratio 2/9 and rises by one order of magnitude from Q 2 = to Q 2 = 100 GeV 2 . 
In Figs. 11 and 12 we compare only pure pomeron contributions to the production cross 
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section; at smaller values of the energy and of the Regge parameter u, the <p / p° and 
(J/^f)/p° production ratios will be further suppressed by the factor f 2 . 

The remarkable restoration of flavor symmetry in the natural scaling variable Q 2 + m\ 
is demonstrated in Figs. 13 and 14, where we present a ratio R(i/k; Q 2 + m v ) of the same 
cross sections taken at equal Q 2 + m v rather than equal Q 2 . A marginal variation of the 
R((J/ty)/p°; Q 2 + my) and R(4>°/p°; Q 2 + m v ) in this scaling variable must be contrasted 
with the variation of the R((J/ty)/ p°; Q 2 = 0) and R((f)°/p°; Q 2 = 0) by the three and one 
orders in the magnitude previously mentioned, respectively, over the same span of Q 2 values 
< Q 2 < 100 GeV 2 . The origin of the slight departures from exact scaling in the variable 
Q 2 + my comes from a well understood difference between the scales At,l and tt,l for 
production of different vector mesons. The same difference of At,l and tt,l brings in the 
energy dependence of R(i/k; Q 2 + m v ). This is a specific prediction from the preasymptotic 
gBFKL dynamics. The radii of the 0° and p° mesons do not differ much and for this reason 
we find a precocious scaling in Q 2 + m v . The energy dependence of the <p°/p ratio also 
turns out very weak. The radii of the p° and J/^f differ much more strongly and the ratio 
R((J/ty)/p°; Q 2 + my) exhibits a somewhat stronger dependence on energy and Q 2 + m v . 
For the same reason, we predict a substantial departure of R(i/k; Q 2 + m v ) from the short- 
distance formula 

which is shown in Fig. 13 by horizontal lines. The formula ( |i~8l) can readily be derived 
generalizing the asymptotic-Q 2 considerations ||, for the further discussion of the crucial 
role of the scaling variable Q 2 + m v in this comparison see below Section 6. 

The case of the uj°,uj' virtual photoproduction is very interesting. Is the p° — u° mass 
degeneracy accidental? Does it imply also similar spatial wave functions in the p and u 
families? The scanning property of diffractive production allows a direct comparison of 
spatial wave functions of the p° and u°. If the u° — p° degeneracy extends also to the spatial 
wave functions, then we predict 

- "°> = i (19) 
a{r -> P°) 9 1 J 

independent of energy and Q 2 . On the other hand, if the radii of the p° and u° are different, 
for instance R w < R p , then the u°/p production ratio must exhibit the Q 2 dependence 
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reminiscent of the </>°/p° ratio. Similarly, a comparison of the u' and p' production can shed 
light on the isospin dependence of interquark forces in vector mesons. 



6 Determination of color dipole cross section from vec- 
tor meson production data 

Inverting Eqs. (0),(0) one can evaluate a(x e ff, r) from the vector meson production data. 
It is convenient to cast Eqs. (|12D , (|T3|) in the form 

2 

mt 



lmM T = g T V4:7ra em C v a(x e ff, r s ) _ 2 " (20) 



VQ 2 m 



2 



lmM L = g L y/ A-Ka em C v a(x eff , r s ) 2 (21) 

In (|20|),(|2lD the coefficient functions gr,i are defined so as to relate the amplitude to a(rs) 
at the well defined scanning radius (fj) with A = 6. The major point of this decomposition is 
that at large Q 2 and/or small rs < Rv, the coefficients qt,l will be very smooth functions of 
Q 2 and energy. The smooth Q 2 and energy dependence of gx,L mostly reflects the smooth and 
well understood Q 2 dependence of scale factors At,l- Such a procedure is somewhat crude 
and the lvnM. T L — cr(x e ff,r s ) relationship is sensitive to the assumed r dependence of the 
dipole cross section a(x e ff,r). Using the dipole cross section || the shape of which changes 
significantly from ui = = 30 up to uj = 3 • 10 6 , we have checked that this sensitivity is 

X eff 

weak. In Fig. 15 we present the Q 2 dependence of the gx,i for different production processes 
at W — 15 GeV and W = 150 GeV. The variation of the resulting coefficient functions 
gx t L from small to large W does not exceed 15%, which is a conservative estimate of the 
theoretical uncertainty of the above procedure. 

The experimentally measured forward cross production section section equals 

Ml *d7 v) ^ = t£ • [ (1 + < t)m2t + e(1 + <^ M "\ ^ 

The difference between ay.i and ay,r for the longitudinal and transverse cross sections and 
the overall effect of the real part is marginal and can safely be neglected compared to other 
uncertainties. Then, making use of the above determined g-r,L and combining Eqs. (|20|) , 



(|21|) and (p2|), we obtain 



1 1 Q 2 + m 2 v 2 / 2 2 Q 2 \ 1/2 



a(x eff , r s ) = - -- — 2 — • lg T + eg L ■ 

J Wly \JOLem \ %, 
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(23) 

t=o 



Here e is the longitudinal polarization of the photon the values of which are taken from 
the corresponding experimental publications. In (^) , (|2"5|) / is the above discussed factor 
which accounts for the non-vacuum Reggeon contribution to the p° production, for <f)° and 
j/vj/ production, / = 1. In the case the experimental data are presented in the form of the 
t- integrated cross section, we evaluate dcr( ^ = Ba tot (j* — > V) using the diffraction 



t=o 



slope I? as cited in the same experimental publication. 



In Fig. 16 we show the results of such an analysis on the low energy [41] and ZEUS [42 



(fp real photoproduction data, on the p° and 0° NMC electroproduction data on the p° 
HERA real and virtual photoproduction (HI f37fl , ZEUS [^, |4(| ^2|), on the fixed target data 
on real photoproduction (EMC |51|, E687 f52"f), on the EMC J/^fr electroproduction data 

0) and on the HERA real photoproduction J/m data (HI gjj, ZEUS |5U|). The error bars 
are the error bars in the measured cross sections as cited in the experimental publications. 
The experimental data on the vector meson production give a solid evidence for a decrease of 
a(x e ff, rs) by one order of magnitude from rs ~ 1.2 fm in 0° real photoproduction down to 
r s « 0.24 fm in the electroproduction of p° at Q 2 = 23 GeV 2 and of J/* at Q 2 = 13 GeV 2 . 
In the region of overlapping values of rs there is a remarkable consistency between the 
dipole size dependence and the absolute values of the dipole cross section determined from 
the data on the p°, 0° and J/\l/ production, in agreement with the flavor independence of 
the dipole cross section. A comparison of determinations of a(x e ff,r) at fixed-target and 
HERA energy confirms the prediction [11], [TJ|, || of faster growth of the dipole cross section 
at smaller dipole size, although the error bars are still large. 

The above determination of a(x e ff,rs) is rather crude for the several reasons. 

1) First, a comparison of the NMC |31| and early EMC data on the p° production 
suggests that the admixture of inelastic process 7*p — > VX could have enhanced the EMC 
cross section by as large a factor as ~ 3 at Q 2 = 17 GeV 2 . The value of <r(7* — > V) thus 
overestimated, leads to cr(x e //, rs) overestimated by the factor ~ a/3, which may be a reason 
why the EMC J/^ electroproduction data |53[j lead consistently to somewhat larger values 
of a(x e ff,rs)- Still, even this factor of ~ a/3 uncertainty is much smaller than the more 
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than the one order of magnitude by which a(x e ff,rs) varies over the considered span of 
rs- In the recent NMC data |31[] a special care has been taken to eliminate an inelastic 
background and the values of a(x e ff,rs) from the p° and 0° production data are consistent 
within the experimental error bars. 

ii) There are further uncertainties with the value of the diffraction slope -8(7* — > V) and 
the curvature of the diffraction cone which affect the extrapolation down to t — 0. The 
experimental situation with the diffraction slopes is quite unsatisfactory; in the case of 
the J/\I/ and of the light vector mesons at large Q 2 , one can not exclude even a ~ 50% 
uncertainty in the value of -8(7* — > V). However, this uncertainty in -8(7* — > V) corresponds 
to < 25% uncertainty in our evaluation of a(x e ff, rs), which is sufficient for the purposes of 
the present exploratory study. 

iii) In addition, there is also the above evaluated conservative < 15% theoretical inaccuracy 
of our procedure. 

iv) Finally, there is a residual uncertainty concerning the wave function of light vector 
mesons. As a matter of fact, if the dipole cross section were known, the diffractive production 
7*p — > Vp would be a unique local probe of the vector meson wave function at r m rs ||; 
this may well become one of the major applications of vector meson production. To this 
aim, the consistency of a(x e ff,r) determined from different reactions indicates that wave 
functions of vector mesons are reasonably constrained by modern spectroscopic models and 
by the leptonic width. 

This is the first direct determination of the dipole cross section from the experimental 
data and our main conclusions on the properties of the dipole cross section are not affected 
by the above cited uncertainties. In Fig. 16 we show also the dipole cross section from the 
gBFKL analysis JH], |(|, which gives a good quantitative description of structure function 
of the photon at small x. We conclude that the color dipole gBFKL dynamics provides a 
unified description of diffractive production of vector mesons and of the proton structure 
function. 

Finally, a comparison of the color dipole analysis of diffractive electroproduction ||, |5], [| 
with the related momentum space analysis of Refs. ||, [7| is in order. At a very large 
Q 2 and/or very short scanning radius, rs <C Ry, the electroproduction probes the wave 
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function of vector mesons and/or the z-distribution amplitude at a vanishing transverse size, 
integrated over the z with the certain ^-dependent factor which emerges in Eqs. (0),©- The 
wave function at the vanishing 3-dimensional separation of the quark and antiquark can be 
related to the width of the leptonic decay, V — *> e + e~ . The form of the z-dependent factor 
is mostly dictated by r| ~ l/(4z(l — z)Q 2 + m 2 / ) which emerges in the integrands of 
after the r integration, and for the asymptotical Q 2 when 4z(l — z)Q 2 ^> my, Brodsky et 
al. [|7| introduced the moment of the longitudinal distribution amplitude 

VV ^dz^ v (r = 0,z) ' 1 j 

One must be careful with the interpretation of r/y , though, because for the very asymmetric 
qq configurations, z(l — z) < M v /Q 2 , the scanning radius stays large even for Q 2 — > oo; 
for instance, precisely these asymmetric configurations dominate the cross section of the 
diffraction dissociation of photons, j*p —>■ Xp, into the continuum states X || [nj. With 
these reservations, we can combine the representations (|20|) , ([2~l|) , the pQCD relationship 
(|TTD and the formula (|]) for the scanning radius, and cast the production amplitude in the 
form (here we focus on the dominant longitudinal amplitude) 



8tt 2 ^ /7 __ _ _ VQt 



M L = — fvV^ctemVvmv— — j—^ ; 2 a s (Q 2 )G(x,Q 2 
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where 



Ifvertvmv ^Q 1 ( Q\ 2 2 

T r s°( r s), (25) 



fv = 7T^nV^e + e-)m v . (26) 



em 



Then, we can present our results for M.^ in terms of this parameter rjy. The first line of 
Eq. ( p5|) gives the asymptotic-Q 2 form of Aii in terms of the gluon structure function of the 
proton, the second line is equivalent to it at large pQCD factorization scale Q 2 and serves as 
a working definition of rjv at moderately large Q 2 and/or moderately small scanning radius 
r$- The above finding that g^ t L only weakly depend on Q 2 and energy, already suggests the 
t]y defined by the second line of Eq. (|2~5"|) will be approximately constant, and now we show 
this is indeed the case. 

Evidently, the resulting values of rjy will depend on the pQCD factorization scale Q 2 . 
The scale parameter r in the pQCD factorization scale Q 2 = t(Q 2 + m v ) was evaluated in 
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0. It is related to the scale parameter B fa 10 in the pQCD formula (|TT| ) and the scale 
parameter A in the scanning radius (0) as r ~ B/A 2 , for production of the longitudinal 
vector mesons in DIS Ref. g finds t l (J/V) « 0.2 and r L (p°) « 0.15. Ref. [0] cites the 
asymptotic leading twist form of (|25|), with neglected in the denominator compared to 



Q 2 and with the pQCD factorization scale Q 2 = Q 2 , besides the more accurate definition 
of the pQCD scale Q 2 we differ from Ref. also by the factor 2 in Eq. flTTl) . The scale A 
in the scanning radius is given by the position of the peak in Wl{Q 2 ,t 2 ), it varies with Q 2 
slightly bringing the slight variation of the scale factors At,l i n > & t large Q 2 it is reasonable 
to take A L (J/^/) = 6 and A L (p°,<j) ) — 8 @. With this choice of A L , our results for the 
production amplitude M. L correspond to the values of rjy shown in Figs. 17 and 18. For 
the nonrelativistic quarkonium, in which z fa |, Eq. ( p4| ) gives rjy fa 2. The T is a good 
approximation to the nonrelativistic quarkonium and we indeed find rjy fa 2. Taking a 
fixed scale A, we neglected the slight variation of r with Q 2 , which propagates into the 
slight variation of rjy with Q 2 . Because the shape of the color dipole cross section varies 
with x e ff, the scale parameter A varies also with x slightly. Taking the x-independent 
A, we cause the slight mismatch of the ^-dependence of the r.h.s. and l.h.s. of Eq. (|24]) , 
which propagates into the weak x e jf dependence of rjy. In Fig. 17 we show rjy for the fixed 
energy W = 150 GeV relevant to the HERA experiments, here the combined Q 2 and x e ff 
dependence of the scale parameter A contribute to the variations of rjy. The issue of the Q 2 
and x dependence of the pQCD factorization scale Q 2 in ([24]) deserves a dedicated analysis, 
here we only wish to focus on the fact that the so determined r]y exhibits a remarkably 
weak variation with Q 2 and x e //. Furthermore, Fig. 18 shows that the x e //-dependence of 
r]y becomes substantially weaker at larger Q 2 . This testifies to an importance of the Q 2 +my 
as a relevant scaling variable, which absorbs major mass corrections to the Q 2 dependence 
of the production amplitude (see also the discussion of the flavor symmetry restoration in 
Section 5). To this end we wish to notice that the expansion 

1 1 ( . 2m 2 v x 



(Q 2 + m 2 v ) 2 Q 4 V 1+ Q 2 "J (2?) 
corresponds to the abnormally large scale 2m v for the higher twist correction to the produc- 
tion amplitude of leading twist. For the light vector mesons, Brodsky et al. cite estimates 
rjy =3-5, our results in Figs. 17 and 18 are very close to these estimates, as it must be 
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expected because the momentum-space technique of Brodsky et al. and our color dipole 
factorization technique are related by the Fourier-Bessel transform. With the present poor 
knowledge of the large dipole distributions in vector mesons and/or the wave functions of 
vector mesons, the variations of rjy in Figs. 17 and 18 and the range of estimates for r]y in 
indicate the range of uncertainty in predictions leptoproduction amplitudes. 

7 Conclusions 

The purpose of this paper has been the phenomenology of diffractive photoproduction and 
electroproduction of ground state (15) and radially excited (25) light vector mesons in 
the framework of the color dipole picture of the QCD pomeron. In this picture, the Q 2 
dependence of production of the 15 vector mesons is controlled by the shrinkage of the 
transverse size of the virtual photon and the small dipole size dependence of the color dipole 
cross section. Taking the same color dipole cross section as used in the previous successful 
prediction of the low x structure function of the proton, we have obtained a good quantitative 
description of the experimental data on diffractive photoproduction and electroproduction of 
15 vector mesons p°,4>° and J/ty. We have presented the first determination of the dipole 
cross section from these data and found a remarkable consistency between the absolute 
value and the dipole size and energy dependence of the dipole cross section extracted from 
the data on different vector mesons. This represents an important cross-check of the color 
dipole picture. The pattern we found for the energy dependence of the dipole cross section 
is consistent with flavor independence and with expectations from the gBFKL dynamics. 
The color dipole picture leads to the restoration of the flavor symmetry and to novel scaling 
relations between the production of different vector mesons when compared at the same 
Q 2 + my. Such relations are borne out by the available data and will be further tested when 
the higher precision data from HERA will become available. Regarding this (Q 2 + rriy)- 
scaling, perhaps still more interesting are the deviations from scaling, which originate from 
a substantial contribution of the large size dipoles even at very large Q 2 s. 

The second class of predictions concerns the rich pattern of an anomalous Q 2 and en- 
ergy dependence of the production of the V'(25) radially excited vector mesons, which 
depends entirely on the quantum mechanical fact that the 25 wave function has a node 
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which makes these anomalies an unavoidable effect. We find a very strong suppression of 
the V'(2S)/V(1S) production ratio in the real photoproduction limit of very small Q 2 . For 
the longitudinally polarized 2S mesons we find a plausible overcompensation scenario leading 
to a sharp dip of the longitudinal cross section a L (2S) at some finite Q 2 = Q 2 n ~ 0.5 GeV 2 . 
The position Q 2 n of this dip depends on the energy and leads to a nonmonotonic energy 
dependence of cfl{2S) at fixed Q 2 . Regarding the experimental choice between the over- 
compensation and undercompensation scenarios in the HERA experiments, the situation 
looks quite favorable because the sign of the p' production amplitude relative to that of the 
p° can be measured directly by the Soding-Pumplin method. At larger Q 2 , the scanning 
radius becomes shorter, and we predict a steep rise of the 2S/1S cross section ratio, typically 
by one order of magnitude on the very short scale Q 2 < 0.5 GeV 2 in agreement with the 
present indications from the E665 data. The flattening of this 2S/1S ratio at large Q 2 is a 
non-negotiable prediction from the color dipole dynamics. Remarkably, the Q 2 dependence 
of the V production offers a unique possibility of distinguishing between 2S radially excited 
and .D-wave vector mesons. 

Finally, in the color dipole framework, a comparison of the Q 2 dependence of the diffrac- 
tive production of the p° and uj° constitutes a direct comparison of the spatial wave functions 
of the two mesons. A comparison of the Q 2 dependence of the u>' and p' production can 
shed light on the isospin dependence of the interquark forces in vector mesons. 
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Appendix. 



Here we present the parameterization of the wave functions of vector mesons in the 
lightcone mixed (r, z) representation. Due to the fact that small size qq configurations 
become important at large Q 2 , one needs to include the short distance hard QCD gluon 
exchange effects so as to make the electromagnetic form factors consistent with the QCD 
predictions. Here we follow a simple procedure suggested in 0, which uses the relativization 



technique of Refs. [55. 56]. We are perfectly aware of the fact the wave functions of light 
vector mesons are still unknown; in the present exploratory study our major concern is to 
have a parameterization which is consistent with the size of vector mesons as suggested by 
the conventional spectroscopic models and has the short distance behavior driven by the 
hard QCD gluon exchange [|56|j . 

Let the Vqq vertex be TqVy^^q where the vertex function T is a function of the lightcone 
invariant variable |55] 



V 2 = \{M 2 -Aml) (28) 

where M is the invariant mass of the qq system 

„ ml + k 2 

M =l^7y (29) 

k and m q are the transverse momentum and quark mass, and z is a fraction of lightcone 
momentum of the meson carried by the quark (0 < z < 1). In the nonrelativistic limit p 
is the 3-momentum of the quark and we have the familiar relationship between the vertex 
function and the momentum space wave function 

*(p 2 ) oc 2 ^ r . (30) 

The hard gluon exchange Coulomb interaction | as ^ , where d is the 3-dimensional 
quark-antiquark separation and as{d) is the running QCD coupling in the coordinate rep- 
resentation, is singular at the origin, d — > 0, but becomes important only at short distances 
d, much smaller than the radius Ry of the vector meson. For this reason, the hard gluon 
Coulomb interaction can be treated perturbatively. Namely, let ty so f t (d) be the wave func- 
tion of the vector meson in the soft, non-singular potential. Solving the Schrodinger equation 
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at small d to the first order in Coulomb interaction, one readily finds the Coulomb-corrected 
wave function of the form 

d 



tf(d) = tf so/t (d) + tf Wt (0)C , exp 
Here a(d) is the "running Bohr radius" equal to 

a(d) 



2Ca(d) 



(31) 



(32) 



8ma s (d) 

where m = m q /2 is the reduced quark mass. The parameter C is controlled by the transition 
between the hard Coulomb and the soft confining interaction; we treat it as a variational 
parameter. (Similar analysis of the correction to the momentum space wave function for 
the short distance Coulomb interaction is reviewed in |^B|). The 3-dimensional Fourier 
transform of the Coulomb-corrected wave function ([H]) reads 



tf(p) = N { (27ri? 2 ) 3/2 exp 



+ c 4 



(33) 



64a 3 (p 2 ) 77 1 
{l + AC 2 a 2 {p 2 )p 2 ) 2 y 

where a(p 2 ) is still given by (^) with the running as(p 2 ) evaluated in the momentum 
representation. 

The relativistic lightcone wave function \l/(z,k) is obtained from \I/(p) by the stan- 
dard substitution of the light cone expression ( p8| , p9|) for the nonrelativistic p 2 in ( |33| ) 
The relativistic wave function thus obtained gives the correct QCD asymptotics 
(Q 2 )/Q 2 of the vector meson form factor, in perfect correspondence to the familiar hard 



QCD mechanism (for the review see |pB |; the more detailed analysis of form factors will be 
presented elsewhere). Then, the lightcone radial wave function is the Fourier transform 

d 2 k 



ir,z 



(2k) 



;^(z, k) exp(ikr) . 



(34) 



With the conventional harmonic oscillator form of ty so f t (d) we obtain the simple analytical 
formula 



<t>is{r, z) 



= tto(lS)|4*(l-2)V27ri2fcexp 



where a(r) is given by Eq. (p2]), (3 = A/B, and 

C 2 a 2 (r)m 2 

\2 i i 



m 2 q R 2 ls 
\z{l - z) 



exp 



2z(l - z)r 2 

R is 





\m 2 q R\ s ] 


exp 


2 



A 



l + ———±-4C 2 a 2 (r)m 2 
z(l — z) q 
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(35) 



(36) 



B 2 



C 2 a 2 (r) 



(37) 



z(l-z) 

For the IS ground state vector mesons we determine the parameters R\ s and C by the 
standard variational procedure using the conventional linear + Coulomb potential models 



| 23j . We check that the resulting wave function are consistent with the experimentally 
measured width of the V — *> e + e~ decay (see Tab. 1). This is one of the major constraints 
because at very large Q 2 and/or r s <C R v , the electroproduction amplitude is controlled 
by the wave function at the vanishing transverse size. For the heavy quarkonia, we check 
that the radii of the IS states are close to the results of more sophisticated solution of 
the Schrodinger equation [23]. The radius of the p° meson given by our wave function 
is consistent with the charge radius of the pion. Still another cross check is provided by 
Ctot{p°N) discussed in Section 3, which comes out very close to the pion-nucleon total cross 
section. 

The node of the radial wave function of the V'{2S) is expected at r n ~ R v far beyond the 
Coulomb region. For this reason, we only modify the soft component of the wave function 
and take the same functional form of the Coulomb correction as for the IS state: 

2z(l - z)r 2 



02s(r, z) 



#o(2£) |4z(l - ^^vri^exp 



™ 2 gR 2 2S 

\z(\-z 



{l -a 

lQa 3 (r 



1 + m n Rnc - r H ^ r 



C 



AB 3 



1(1 L2S Az(l-z) 
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exp 
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~™ 2 q Rls 


exp 
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(38) 



The new parameter a controls the position r n of the node. The two parameters a and R2S 
are determined from the orthogonality condition 



Nr. 



dz 



2vrJo z 2 (l 



d 2 r 



[mfyiir, z)(f) k (r, z) + [z 2 + (1 - z^ldr&ir, z)][d r (f) k (r, z)}} = 5 ik 



(39) 



and from the 25 — 15' mass splitting evaluated with the same linear+ Coulomb potential. For 
the heavy quarkonia, we can check the resulting V'(2S) wave function against the accurate 
data on the width of the V'{2S) — > e + e~ decay, the agreement in all the cases is good. The so 
determined parameters, the quark masses used and some comparisons with the experiment 
are summarized in Table 1. It is Ca(r) which defines at which radii the interaction is 
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important. A posteriori, for light vector mesons C is found small, the radius Ca(r) is indeed 
small and the resulting parameters are consistent with the assumption that the Coulomb 
interaction is a short-distance perturbation. Furthermore, for the light vector mesons we 
find Ris ~ i?25- For heavier mesons C is larger and Coulomb effects are becoming more 
important and i?2S > Ris m the ratio closer to the one for the Coulomb system (see Table 1). 
Because our Ansatz for the relativistic wave function has the correct short-distance QCD 
behavior and gives a reasonable description of soft cross sections, we believe it provides 
a reasonable interpolation between the soft and hard regimes in the electroproduction of 
vector mesons. 



31 



References 

[1] A.Donnachie and RV.LandshofT, Phys. Lett. B185 (1987) 403; J.R.Cuddell, Nucl. 
Phys. B336 (1990) 1. 

[2] B.Z.Kopeliovich and B.G.Zakharov, Phys. Rev. D44 (1991) 3466. 

[3] M.G.Ryskin, Z. Phys. C57 (1993) 89. 

[4] B.Z.Kopeliovich, J.Nemchik, N.N.Nikolaev and B.G.Zakharov, Phys. Lett. B309 

(1993) 179. 

[5] B.Z.Kopeliovich, J.Nemchik, N.N.Nikolaev and B.G.Zakharov, Phys.Lett. B324 (1994) 
469. 

[6] J.Nemchik, N.N.Nikolaev and B.G.Zakharov, Phys.Lett. B341 (1994) 228. 
[7] S.J.Brodsky et al., Phys. Rev. D50 (1994) 3134. 
[8] J.R.Forshaw and M.G.Ryskin, Z. Phys. C (1995) to be published 
[9] N.N. Nikolaev and B.G. Zakharov, Z. Phys. C49 (1991) 607; Z. Phys. C53 (1992) 331; 
[10] N.Nikolaev and B.G.Zakharov, JETP 78 (1994) 598; Z. Phys. C64 (1994)631; 

[11] N.N.Nikolaev, B.G.Zakharov and V.R.Zoller, JETP Letters 59 (1994) 8; JETP 78 

(1994) 866; Phys. Lett. B328 (1994) 486. 

[12] A.E.Mueller, Nucl. Phys. B415 (1994) 373; A.E.Mueller and H.Patel, Nucl. Phys. 
B425 (1994) 471; Z.Chen and A.E.Mueller, Nucl. Phys. B451 (1995) 579. 

[13] N.N.Nikolaev, Comments on Nucl. Part. Phys. 21 (1992) 41. 

[14] N.N.Nikolaev and B.G.Zakharov, Phys. Lett. B327 (1994) 157. 

[15] J.Nemchik, N.N.Nikolaev and B.G.Zakharov, Phys.Lett. B339 (1994) 194. 

[16] NMC Collab., P.Amaudruz, et al., Nuc. Phys. B 371 (1992) 3. 



32 



[17] P.Frabetti et al., presented on Int. Europhys. Conf. on HEP, Brussels, July 27 - August 
2, 1995. 

[18] PSoding, Phys. Lett. 19 (1966) 702. 

[19] J.Pumplin, Phys. Rev. D2 (1970) 1859. 

[20] K.Abe et al., Phys. Rev. Lett. 53 (1984) 751. 

[21] N.N.Nikolaev and B.G.Zakharov, Phys. Lett. B327 (1994) 149. 

[22] F.E.Low, Phys. Rev. D12 (19750 163; J.F.Gunion and D.E.Soper, Phys. Rev. D.15 
(1977) 2617; A.B.Zamolodchikov, B.Z.Kopeliovich and L.I.Lapidus, JETP Letters 33 
(1981) 612; G.Bertsch, S.Brodsky, A.S.Goldhaber and J.F.Gunion, Phys. rev. Lett. 47 
(1981) 297. 

[23] W.Buchmuller, S.-H.H.Tye, Phys. Rev. D24 (1981) 132 
[24] J.D.Bjorken, J.B.Kogut and D.E.Soper, Phys. Rev. D3 (1971) 1382. 
[25] S.R.Gevorkyan and A.V.Tarasov, Dubna preprint E2-94-465 (1994). 
[26] V.N.Gribov, A.A.Migdal, Sov. J. Nucl. Phys. 8 (1969) 703 

[27] V.Barone, M.Genovese, N.N.Nikolaev, E.Predazzi and B.G.Zakharov, Z. Phys. C58 
(1993) 541; Int.J.Mod.Phys A, 8 (1993) 2779. 

[28] N.N.Nikolaev and B.G.Zakharov, Phys. Lett. B332 (1994) 184. 

[29] T.H.Bauer, R.D.Spital, D.R.Yennie and F.M.Pipkin, Rev. Mod. Phys. 50 (1978) 261. 

[30] E665 Collab., M.R.Adams et al, Phys. Rev. Lett. 74 (1995) 1525 

[31] NMC Collab., M.Arneodo, et al, Nucl.Phys. B 429 (1994) 503 

[32] ZEUS Collab., M.Derrick et al., Phys. lett. B356 (1995) 601 

[33] S.J.Brodsky and A.H.Mueller, Phys. Lett. B206 (1988) 685. 

[34] M.R.Adams et al., Phys. Rev. Lett. 74 (1995) 1525. 

33 



[35] H.Schellman, Blois Workshop on Diffractive and Elastic Scattering, Blois, June 1995 

[36] A.Donnachie and P.V.L.Landshoff, Phys. Lett. B296 (1992) 227; Phys. Lett. B348 
(1995) 213 

[37] HI Collab., T.Ahmed et al., Exclusive p° Production in Deep-inelastic Scattering Events 
at HERA, presented on Int. Europhys. Conf. on HEP, Brussels, July 27 - August 2, 
1995, paper EPS-0490 

[38] OMEGA Collab., D.Aston et al. Nucl. Phys. B209 (1982) 56; J.Park et al. Nucl. Phys. 
B36 (1972) 404; R.M.Egloff et al. Phys. Rev. Lett. 43 (1979) 657 

[39] ZEUS Collab., M.Derrick et al., Z. Phys. C63 (1994) 391. 

[40] ZEUS Collab., M.Derrick et al., DESY 95-143, submitted to Z. Phys. C (1995) 

[41] J. Busenitz et al., Phys. Rev. D40 (1989) 40 and refrences therein. 

[42] ZEUS Collab., M.Derrick et al., Elastic Photoproduction of u, <fi and p' mesons at 
HERA, presented on Int. Europhys. Conf. on HEP, Brussels, July 27- August 2, 1995, 
paper EPS-0389 

[43] O.Benhar, B.G.Zakharov, N.N.Nikolaev and S.Fantoni, Phys. Rev. Lett. 74 (1995) 3565 

[44] A.Donnachie and A.B.Clegg, Z. Phys. C42 (1989) 663 and references therein. 

[45] Review of Particle Properties, Phys. Rev. D50 (1994) 1177 

[46] E665 Collab., C.Salgado, Workshop on CT at CEBAF, May, 1995 

[47] N.N.Nikolaev, B.G.zakharov and V.R.Zoller, Phys. Lett. B366 )1996) 337. 

[48] J.Nemchik, N.N.Nikolaev, E.Predazzi and B.G.Zakharov, BFKL pomeron and diffrac- 
tive electroproduction of heavy vector mesons at HERA, paper in preparation 

[49] HI Collab., T.Ahmed et al., Phys. Lett. B338 (1994) 507 

[50] ZEUS Collab., M.Derrick et al, Phys. Lett. B350 (1995) 120 



34 



[51] EMC Collab., J.J.Aubert et al, Phys. Lett. B89 (1980) 267. 

[52] L.Frabetti et al, Phys. Lett. B316 (1993) 197. 

[53] EMC Collab., J.J.Aubert et al, Nucl. Phys. B213 (1983) 1; 
J.Ashman et al.,Z. Phys. C39 (1988) 169. 

[54] EMC Collab., J.J.Aubert et al, Phys. Lett. B161 (1985) 203. 

[55] M.V.Terentiev, Sov. J. Nucl. Phys. 24 (1976) 106. 

[56] S.J.Brodsky and G.P.Lepage, in Perturbative Quantum Chromo dynamics, editor 
A.H.Mueller, World Scientific, Singapore, 1989. 



35 



Figure captions: 

Fig. 1 - The color dipole model predictions for the total cross section a to t{VN) for the 
interaction of the light vector mesons p, p',<f) and 0' with the nucleon target as a 
function of c.m.s. energy W. 

Fig. 2 - The color dipole model predictions for the Q 2 dependence of the observed cross 
section a (7* — > V) = o"t(7* — > V) + £0^(7* — > V) of exclusive p° and 0° production 



vs. the low-energy NMC [31] and high-energy ZEUS [32] and HI |37| data. The top 
curve is a prediction for the p° production at W = 70 GeV, the lower curves are for 
the p°, 0° production at W = 15 GeV. The dashed curve for the p° shows the pure 
pomeron contribution 0^(7* — > p°), the solid curve for the p° shows the effect of 
correcting for the non-vacuum Reggeon exchange as described in the text. 

Fig. 3 - The color dipole model energy dependence predictions for forward real photoproduc- 
tion of p° mesons compared with fixed target data |38| and high energy datum from the 



ZEUS experiment at HERA collider |39], |40|| . The dashed curve is the pure pomeron 



exchange contribution, the solid curve shows the correction for the the non-vacuum 
Reggeon exchange as described in the text. 

Fig. 4 - The color dipole model predictions for the energy dependence of real photoproduc- 
tion of the 0° mesons compared with fixed target [[11]] and high energy ZEUS data 
(open square for the 0° [Q, solid circle for the p° [^. [l(| ). 

Fig. 5 - The color dipole model predictions of the forward differential cross sections da^^il* — 
V)/dt\ t =Q for transversely (T) (top boxes) and longitudinally (L) (middle boxes) polar- 
ized p° and 0° and for the polarization-unseparated da^* — > V)/dt\ t=Q = daxi'y* — > 
V)/dt\ t =o + ed(Ti(7* — > V)/dt\ t =o (bottom boxes) for e = 1 as a function of the c.m.s. 
energy W at different values of Q 2 . 

Fig. 6 - The color dipole model predictions for the Q 2 and W dependence of the ratio of 
the longitudinal and transverse differential cross sections in the form of the quantity 
Rlt = ^T dCTy^*^!/) , where my is the mass of the vector meson. The solid and dashed 
curves are for W = 15 GeV and W = 150 GeV. 
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Fig. 7 - The color dipole model predictions for the dependence on the scaling variable Q 2 + 
uiy of the polarization-unseparated da(-y* — > V)/dt\ t=0 = ^0^(7* — > K)/<it| t=0 + 
eda L (Y V)/dt\ t=0 for e = 1 at the HERA energy W = 100 GeV. 

Fig. 8 - The color dipole model predictions for the Q 2 and W dependence of the ratios 
a(Y -> p'(2S))/a(Y -> p°) and a( 7 * -> (f/(2S))/a(Y -> 0°) for the (T) and (L) 
polarization of the vector mesons. 

Fig. 9 - The color dipole model predictions for the Q 2 dependence of the ratio of the 
polarization-unseparated forward production cross sections <icx(7* — > p'(2S))/da( , -f* —>■ 
p°) and g?<t(7* — > (/>' (2S))/ dafa* — > 0°) for the polarization of the virtual photon e = 1 
at the HERA energy W = 100 GeV. 

Fig. 10 - The color dipole model predictions of the forward differential cross sections da^^il* ~ 
V')/dt\t=o for transversely(T) (top boxes) and longitudinally (L) (middle boxes) po- 
larized radially excited vector mesons p'{2S) and <p'(2S) and for the polarization- 
unseparated da(j* — > V')/dt\ t =o = da T (j* — > V')/dt\ t=0 + eda L (j* — > V')/dt\ t=0 for 
e = 1 (bottom boxes) as a function of the c.m.s. energy W at different values of Q 2 . 

Fig. 11 - The color dipole model predictions for the energy dependence of the ratio of the 
polarization-unseparated forward production cross sections dcr('-f* — > (p )/da('~f* — > p°) 
for the polarization of the virtual photon e = 1 at different values of Q 2 . 

Fig. 12 - The color dipole model predictions for the energy dependence of the ratio of the 
polarization-unseparated forward production cross sections da{j* — > J/^)/da(j* — > 
p°) for the polarization of the virtual photon e = 1 at different values of Q 2 . 

Fig. 13 - Approximate scaling in the variable Q 2 + rriy for the ratio of the polarization- 
unseparated forward production cross sections da{pf* — > 0°)/dcr(7* — > p°) andrfcr(7* — ■> 
J/^f) jdoi^f — > p°) for the polarization of the virtual photon e — 1. The horizontal 
dotted straight lines show the ratio corresponding to Eq. (|18|). 

Fig. 14 - Approximate scaling in the variable Q 2 + my for the ratio of the polarization- 
unseparated forward production cross sectionsf da (7* — > (p )/da{j* — > p°) and<icr(7* — ; 
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J/^)/da(j* — > p°). at c.m.s. energy W = 150 GeV (the polarization of the virtual 
photon e = 1). 

Fig. 15 - The Q 2 dependence of the coefficient functions g^L at W — 15 GeV (dashed curve) 
and W = 150 GeV (solid curve). 

Fig. 16 - The dipole size dependence of the dipole cross section extracted from the experimen- 
tal data on photoproduction and electroproduction of vector mesons: the NMC data 
on (f)° and p° production f3~I| , the EMC data on J/ty production |53|, the E687 



data on J/ty production f5"2| , the FNAL data on p° production ||41|| , the ZEUS data 
on 4>° production f42|, the ZEUS data on p° production |39| f|(| [3^1 , the HI data on p° 



production |37j and the average of the HI and ZEUS data on J/^f production |49| , ^0[ 



The dashed and solid curve show the dipole cross section of the model |^T[ evalu- 
ated for the c.m.s. energy W = 15 and W = 70 GeV respectively. The data points at 
HERA energies and the corresponding solid curve are multiplied by the factor 1.5. 

Fig. 17 -The Q 2 dependence of the parameter r]y in the representation (25) for the amplitude 
of leptoproduction of different vector mesons at fixed energy W = 150 GeV. 

Fig. 18 -The x e ff dependence of the parameter r]y in the representation (|25|) for the amplitude 
of leptoproduction of different vector mesons at several values of Q 2 . 
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parameter 


P° 


P' 


<f)° 


4>' 


J/tf 




T 




R 2 [fm 2 ] 


1.37 


1.39 


0.69 


0.83 


0.135 


0.248 


0.015 


0.047 


C 


0.36 


0.28 


0.53 


0.44 


1.13 


0.99 


1.99 


1.72 


m q [GeV] 


0.15 


0.15 


0.30 


0.30 


1.30 


1.30 


5.00 


5.00 


a 




0.86 




0.94 




1.20 




1.53 


Rv [fm] 


1.30 


2.28 


0.91 


1.68 


0.41 


0.83 


0.19 


0.42 


Am(2S - IS) [GeV] 




0.73 




0.64 




0.60 




0.55 


r(e + e~) [keV] 


6.29 


2.62 


1.23 


0.47 


5.06 


1.78 


1.20 


0.54 


T ex P(e + e-) [keV] 


6.77 




1.37 


0.48 


5.36 


2.14 


1.34 


0.56 




±0.32 




±0.05 


±0.14 


±0.29 


±0.21 


±0.04 


±0.140 



Table 1: The parameters R 2 , C, m q and a of the vector mesons wave function and some 
of the observables evaluated with these wave functions: the r.m.s. Ry, the leptonic width 
r(e + e~) and the V'(2S) — V(1S) mass splitting. The values of T(e + e~) from the Particle 
Data Tables p5|/ is shown for the comparison. 
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